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SUMMARY 

This  paper  reports  on  calculations  of  the  electron 
energy  distribution  occurring  in  gases  and  gas 
mixtures  used  in  switching.  For  externally 
ionized  high  pressure  discharges  the  optimization 
of  the  conductivity  as  a  function  of  a  gas  mixture 
(usually  a  Rarnsauer  gas  plus  a  molecular  gas)  and 
of  E/N  is  examined.  The  power  loading  and  pulse 
length  change  the  vibrational  populations  of  the 
molecular  gas  and  may  reduce  the  conductivity 
unless  the  vibrational  relaxation  is  sufficiently 
rapid.  For  high  repetition  rates  we  also  examine 
the  influence  of  dissociation  on  the  plasma 
conductivity  and  how  the  decay  processes  might  be 
influenced  by  the  degree  of  vibrational 
excitation. 

At  low  pressures  it  is  more  appropriate  to  treat 
the  switch  plasma  as  a  negative  glow  plasma.  At 
low  electron  densities,  this  requires 
consideration  of  the  ionization  and  excitation  by 
the  fast  electrons  from  the  cathode  sheath.  At 
high  electron  densities  the  electron  distribution 
function  tends  to  a  Maxwellian  however  some 
electrons  with  energies  greater  than  the  injection 
energy  appear  due  to  velocity  diffusion.  These 
features  must  be  included  in  the  traditional 
treatment  of  thyratrons  and  other  switches. 


Since  the  electron-electron  Coulomb  collision 
frequency  varies  roughly  as  n  v  ,  these 
collisions  are  especially  effective  in  the  elastic 
regime  where  there  are  no  severe  energy  loss 
processes,  and  they  may  have  an  influence  at 
fractional  ionization  as  low  as  10-/  for  noble 
gases. 

The  switch  plasmas  fall  into  two  main  categories: 
high  pressure  and  low  pressure.  The  approaches 
for  calculating  the  electron  kinetics  of  these 
discharges  are  significantly  different.  At  high 
pressures  the  plasma  can  still  be  diffuse.  This 
is  made  possible  by  external  electron  beam 
ionization.  Otherwise  an  arc  will  form  or  indeed 
it  may  still  develop  if  the  pulse  length  is  too 
long.  Criteria  for  this  to  occur  have  been 
developed  using  streamer  theory  or  ionization  - 
thermal  instability  analyses5.  In  the  fully 
developed  arc,  equilibrium  plasma  theory  can  be 
applied0.  The  spark  gap  is  described  by  this 
approach  and  its  principal  limitation  of  low 
repetition  rate  is  determined  by  the  collisional 
coupling  term  between  the  atoms  and  the  electrons. 
This  coupling  maintains  conductivity  even  when  the 
electric  field  is  removed.  This  means  that  the 
electron  density  decays  slowly  unless  the  gas 
travels  across  the  discharge  gap  very  rapidly. 


INTRODUCTION 


In  order  to  calculate  the  electron  conductivity  of 
a  plasma  switch  it  is  necessary  to  measure  or  to 
be  able  to  calculate  the  electron  spatial  and 
velocity  distribution  f  (y,  r,t).  The  electron 
density  and  current  are  then 


(2) 

lik 


The  electron  drift  velocity  W„  =L3j'l&.  This 
velocity  is  generally  small  compared  to  the 
electron  random  velocity  v.  The  electron  energy 
distribution  function  (EEDF)  f  (v,r,t)  is  de¬ 
scribed  by  the  Boltzmann  transport  equation1. 
This  continuity  equation  in  phase  space  has  the 


form 


if 


where  the  subscripts  on  the  divergence  operators 
denote  the  independent  variable;  v^  is  the  particle 
velocity  and  ja  is  the  particle  acceleration.  The 
collision  term  includes  the  effects  of  elastic 
collisions,  inelastic  collisions  with  atoms  or 
molecules  and  electron-electron  collisions.  The 
collisions  are  usually  treated  as  source  terms, 
instantaneously  transferring  electrons  from  one 
volume  of  velocity  space  to  another.  The  elec¬ 
tron-electron  collisions  will  produce  a  Maxwellian 
velocity  distribution  when  the  fractional 
ionization  is  sufficientl^-large. ,  Ihis  is 
generally  greater  than  10  3  -  10-H  ’  . 


ELECTRON  KINETICS  IN  THE 
HIGH  PRESSURE  DIFFUSE  DISCHARGE 

In  high  pressure  diffuse  discharges  that  may  be 
suitable  for  on/off  switching  the  EEDF  is  usually 
dominated  by  electron  neutral  collisions.  The 
collision  term  is  then  given  by  the  Boltzmann 
collision  integral  and  a  knowledge  of  the  larger 
elastic  and  inelastic  collision  processes  is 
required.  If  the  applied  field  is  not  too  large 
(i.e.  the  energy  gained  between  collisions  is 
small  compared  to  the  mean  electron  energy), 
perturbation  methods  for  solving  the  Boltzmann 
transport  equation  may  be  applied.  The 
distribution  function  is  normally  expanded  in 
spherical  harmonics  in  velocity  space.  The  first 
term  represents  the  isotropic  distribution  and  the 
higher  order  terms  represent  the  anisotropy 
induced  by  the  perturbing  fields. 


P.  (0)  is  ith  Legendre  polynomial 

0  is  polar  angle  w.r.t.  electric  field 
then  (3) 
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The  Boltzmann  equation  may  be  rewritten  in  the 
form  of  a  system  of  coupled  equations 


Where  ^m  is  the  momentum  loss  collision  frequency 


?  =  N  aj- 

An  o 


integrated  over  all  solid  angles  andu/f4-,  (7/is  the 
differential  cross-section  for  electron  scattering. 


COS&)A3Q  (7) 

iQkoh 


If  the  applied  electric  field  is  not  changing 


Note  that  in  the  absence  of  electron  density 
gradients , 

Pfo  (9) 


/  ^9 
yyyx 


This  remark  is  made  here  so  as  to  interpret 
calculations  discussed  below.  If  the  principal 
cross-sections  are  available  then  it  is  possible 
to  derive  the  breakdown  condition,  the  drift 
velocity  and  other  transport  properties.'  There 
is  some  improvement  fog  high  fields  by  going  to  a 
MontegCarlo  simulation0,  or  using  higher  order 
terms  ,  however  the  spatially  averaged  Townsend 
coefficient  is  a  reasonable  approximation  except 
in  the  cathode  fall. 


All  of  the  above  approaches  have  been  cultivated 
very  carefully  and  used  very  successfully  by  those 
interested  in  deriving  data  from  swarm 
experiments.  However  for  discharges  of  switching 
interest  there  are  other  effects  that  must  be 
considered  and  that  cause  significant  correction 
to  the  plasma  conductivity.  First  of  all  the 
switch  will  have  a  much  higher  electron  density 
than  the  drift  tube.  Thus  electron-electron 
colfisions  must  be  considered.  Next,  in  practical 
situations  the  gas  composition  is  not  independent 
of  the  excitation  conditions.  Therefore  either 
deliberately  or  inherently  because  of  dissociation 
one  has  a  gas  mixture.  These  features  lead  to 
some  new  effects,  some  of  which  might  be 
exploitable  in  pulse  power  switching.  The  third 
difference  is  that  in  the  discharge  any  diatomic 
gas  will  be  vibrationally  excited  and,  especially 
if  the  gas  is  homonuclear,  the  excited  vibrational 
levels'  populations  can  become  a  significant 
percentage  of  the  total  gas  density.  The 
vibrational  populations  have  two  influences: 
first  they  can  change  many  of  the  electron  kinetic 
rates.  Some  of  the  excitation  or  ionization 
thresholds  are  altered  (usually  lowered)  and  in 
some  cases  the  rate  constant  is  dramatically 
modified  (e.g.  dissociative  attachment  in  IL , 

D,,  and  HCL):  The  second  influence  is  that  a 
vTbrational  population  reacts  back  on  the  electron 


gas  through  collisions  of  the  second  kind.  This 
means  the  electron  energy  distribution  at  a  given 
E/N  for  a  power  loaded  gas  is  quite  different 
from  that  for  present  drift  tube  experiments.  (In 
principle  some  drift  tube  experiments  in 
vibrationally  excited  gases  are  possible,  but  no 
one  has  attempted  to  exploit  the  long  lifetime  of 
diatomic  molecules  in  this  way). 

The  Boltzmann  transport  equation  calculations  are 
in  reasonable  order  for  the  pure  rare  gases, 
hydrogen,  oxygen  and  several  laser  gases.  Data 
for  mercury  is, also  available  . An  example  for 
oxygen  plasmasiU  is  shown  in  fig  1.  Note  that 
these  distirbutions  cannot  be  approximated  by 
either  a  Maxwellian  or  a  Druyvesteyn 
distribution.  In  hydrogen  however,  a  Maxwellian 
distribution  is  a  reasonable  approximation  at  low 
electron  densities  and  low  power  loading. 


However  these  analyses  assume  that  the  discharge 
gases  are  not  altered  by  the  power  loading. 

Under  high  power  switch  conditions  this 
assumption  often  breaks  down.  If  the  original 
gas  is  molecular, .then  it  will  dissociate  with  a 
time  constant  where 

At  high  powers  this  effect  will  tend  to 
compensate  the  gas  depletion  caused  by  ionization 
however  it  is  necessary  to  allow  for  the  fact 
that  the  species  has  changed.  This  becomes  very 
important  for  the  slower  electrons  as  the 
effective  vibrational  cross-sections  in  a 
molecular  gas  are  being  reduced  as  dissociation 
proceeds.  There  are  no  equivalent  loss  processes 
in  the  same  energy  range  until  the  excited  state 
densities  become  very  high.  Our  calculations  for 
hydrogen  indicate  that  the  dissociation  rate 
exceeds  the  ionization  rate  by  approximately  an 
order  of  magnitude  over  most  of  the  E/N  range. 
Therefore  the  dissociation  will  occur  faster  than 
the  ionization  growth  and  must  be  considered  in 
most  pulsed  discharges. 

Some  examples  of  the  effects  of  dissociation  on 
the  EEDF  are  shown  in  fig  2  (oxygen)  and  in  fig  3 
(hydrogen)  Mote  that  the  rates  for  processes 
controlled  by  the  fast  electrons  will  be 
significantly  altered.  The  excitation  rates 
derived  from  the  EEDF's  for  molecular  hydrogen 
are  shown  in  fig  4.  At  E/N  values  greater  than 
20  Td  a  large  fraction  of  the  discharge  energy  in 
hydrogen  or  in  deuterium  goes  into  vibrational 
excitation.  This  then  can  change. the  excitation, 
ionization  and  dissociation  rates11  and  the 
EEDF  is  altered  by  the  effects  mentioned  above. 
Under  some  conditions  it  is  believed  the  process 
can  be  runaway  and  lead  to  instability  even 
earlier  than  thermal  heating.  Most  dramatic,  the 
dissociative  attachment  rate  in  hydrogen  and 
deuterium  is  predicted  to  be  changecLorders  of 
magnitude  on  vibrational  excitation  .  The 
calculated  results  for  hydrogen  are  shown  in 
figure  5. 

Generally,  electron  drift  velocity  experiments 
have  concentrated  on,simple  gases.  However 
Cottrell  and  Walker1,3  in  1965  did  measure  the 
electron  drift  velocity  in  a  large  number  of 
polyatomic  gases.  Some  of  these  gases  showed 
remarkably  high  drift  velocities  at  low  E/N 
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values  -  just  the  ideal  conditions  for  a  switch  ELECTRON  KINETICS  IN  THE 

in  the  on-condition.  LOW  PRESSURE  DISCHARGE 


For  low  energy  electrons,  if  we  assume  isotropic 
scattering,  then  on  average,  the  electron  motion 
in  the  field  direction  is  reduced  to  zero  in  a 
collision.  Then  K4,  =  £ea  [  2smAT 


or 


% 


-AT 


(10) 


where  X  =  electron  mean  free  path  =  i/A/Qj- 

Ee.  l2/Yn,NC<-y- 

In  steady  estate  'U /k  =  a  Vy^" 
where  is  the  average  fractional  energy 

loss  in  a  collision  and  a  is  a  constant  of  the 
order  unity  then  7  .  U  )JL 


H, 
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(11) 


While  arid  Qj( total  cross-section)  are  func¬ 

tions  of  E/N,  it  is  apparent  that  a  high  drift 
velocity  is  favored  by  a  large  fractional  energy 
loss  and  a  small  total  cross-section  (Compare 
eqn.  (9) 


by  a 
loss 


1  arge 
and  a 


a  large  drift  velocity  is  favored 
AJ~  (caused  by  a  large  energy 
Q/yn/,  electron  momentum 
transfer  cross-section).  Thus  if  a  gas  has  a 
Ramsauer  minimum  it  satisfies  one  of  these 
conditions;  normally  the  heavier  rare  gases  do  not 
have  inelastic  losses  near  the  energy  of  the 
Ramsauer  minimum  and  so  do  not  have  exceptionally 
hjgh  electron  drift  velocities.  Gases  lik eCRn, 
C2r4’and  appear  to  be  Ramsauer  gases  and^ 

also  have  low  energy  vibrational  modes  that  are 
excited  by  electronic  collisions,  and  so  would 
qualify.  It  is  possible  to  synthesize  a  gas  with 
desirable  high  mobility  features  by  using  a  rare 
gas-molecular  mixture.  Some  of  our  calculated 
results  are  shown  in  figure  6  for  krypton  -  carbon 
monoxide  mixtures.  Experimental  measurements 
have  recently  confirmed  these  models .for  nitrogen 
-  argon  and  silane  -  argon  mixtures. 


In  an  actual  discharge,  if  the  vibrational 
manifold  becomes  saturated,  the  net  energy  loss  is 
reduced  as  there  is  energy  transfer  back  from  the 
excited  states.  This  will  reduce  the  enhancement 
of  the  electron  drift  velocity.  Therefore  a_ 
triple  mixture  with  a  species  to  relax  the 
vibrational  states  may  be  necessary  if  the  energy 
loading  itself  does  not  produce  sufficient 
dissociation  to  cause  the  relaxation.  (This  may 
partially  account  for  some  of  the  different 
conductivities  reported  for  Ar  N2  mixtures  by 
different  experiments). 


For  switches  such  as  thyratrons  that  operate  on 
the  low  pressure  side  of  the  Paschen  curve,  the 
derivation  of  the  electron  kinetics  is  more 
complicated  than  in  the  high  pressure  case.  Many 
of  these  discharges  are  better  characterized  as 
negative  glows  rather  than  positive  column 
discharges.  Thus  E/N  is  not  as  unique  a  scaling 
parameter.  In  many  of  the  discharges  the  EEDF  is 
not  in  equilibrium  with  the  local  value  of  the 
electric  field.  Long15  and  Segur  et  al. 
have  treated  the  non-equilibrium  regions  of  a 
weakly  ionized  gas.  These  analyses  are  very 
useful  in  describing  the  low-pressure  cold  cathode 
discharge.  (Actually  the  non-equilibrium  also 
will  occur  in  the  high  pressure  discharge,  however 
the  spatial  extent  is  extremely  small). 

Conscious  of  the  important  of  non-equilibrium 
effects,  let  us  examine  the  hot  cathode  discharge 
which  appears  to  be  more  controllable  as  a  switch 
and  has  attracted  renewed  interest.  A  Boltzmann 
type  analysis  can  be  used  for  the  case  where  the 
mean  free  path  is  no  longer  small  compared  to  the 
discharge  characteristic  lengths.  A  very  useful 
approach  to  this  discharge  was  introduced  bv 
K.  R.  Emeleus  and  applied  by  Leckey  et.al.1  . 
Additional  studies  have  been  made  by 
Aleskovskii1  .  The  cathode  sheath  injection  of 
fast  electrons  plays  an  important  role  in  the  the 
thyratron  discharge  where  the  pessure  is  about  0.3 
torr  and  the  mean  free  path  is  about  C.3cm.  Thus 
the  distribution  function  of  ionizing  electrons  is 
given  as  a  fixed  condition  along  the  cathode 
boundary  of  the  negative  glow  (NG).  In  this  case 
the  diffusion  of  fast  electrons  in  the  NG  can  be 
described  in  the  P..  approximation  by  the  scalar 
flux  equation  _ 

I  (A  E)  =  E/  Q  ( i2) 

where  Sis  the  unit  velocity  vector,  E  is  the 
electron  energy,  £  is  the  position  vector  and  f 
(r,  E,i4)  is  the  distribution  function.  The 
inelastic  collisions  are  loss  terms  and  the 
diffusion  of  electrons  from  the  known  boundary 
source  is  the  gain  term.  The  kinetic  equation  is 

-  Q.Mlk.s)  os) 

+2.  Q:« (ei-E* )l(ij €+  £  *)  -  o 

where  De  is  the  free  electron  diffusion  coefficient 
E.  -  excitation  energy 


The  EEDF  in  CO-Krypton  mixtures  is  shown  in  Fig  7. 
This  gas  mixture  demonstrates  a  remarkable  drift 
velocity  enhancement  over  either  pure  gas.  The 
velocity  space  diffusion  control  of  IL  is  shown 
by  the  curves  at  low  E/N  where  the  steep  change  in 
the  isotropic  distribution  with  energy  is  caused 
by  the  inelastic  collisions  with  CO  at  energies 
near  the  krypton  momentum  transfer  minimum.  Thus 
there  is  correspondence  between  the  numerical 
results  and  the  analytical  expression  (eqn  9). 


The  condition  that  the  cathode  sheath  voltage  is 
not  much  greater  than  the  ionization  potential  is 
usually  satisfied  in  a  hot  cathode  discharge. 

then  D.eV%-e(fast)--  X /7e  (^t)=  O  (14) 

where  K  =  NQ.^1/- 


The  general  solution  is  expressed  by  a  linear 
summation  of  the  form1 

V?^5) 


n.jAy)^[c,vexp(?vf)  +C2^ 

*  \Jo(evA) 

Where  V  20  was  needed  to  describe  the 
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immediate  cathode  region.  The  contours  of  the 
fast  electron  density  illustrated  in  Fig  8  are 
obtained  from: 


The  spatial  distributions  of  the  ionization  and 
excitation  rates  can  then  be  calculated.  The 
remainder  of  the  discharge  potential  is 
distributed  across  the  glow  and  its  spatial 
distribution  is  determined  by  the  discharge 
current  continuity  requirement.  The  low  energy 
portion  of  the  EEDF  is  assumed  to  be  in 
equilibrium  with  these  E/N  values. 


At  higher  current  densities  the  fast  electron 
energy  loss  will  be  dominated  by  Coulomb 
collisions.  The  Boltzmann  equation  then  reduces 
to  the  Fokker  -  Planck  equation.  Hiskes  and 
KaroiS  obtain 


(17) 

where  f  is  the  EEDF  for  the  slowing  electrons,  q 
is  the  electron  drag  rate  and  J  is  the  injected 
current/unit  energy.  The  electron  drag  is 
evaluated  assuming  that  it  is  dominated  by  the 
near  Maxwellian  slow  group.  If  the  injected 
current  is  constant  over  a  range  of  energy 


E  <~  E p  < 


—(Ef  E0)T 

-e  ~%T~}\ 


-(E-Ef) 
e  ~TT 


E  >Ef 


Hiskes  and  Karo  have  applied  this  analysis  to  the 
Berkeley  56kw  hydrogen  ion  source.  Fig  9  shows 
the  calculated  distribution  for  this  low  pressure 
hydrogen  discharge. 


Arguments  have  usually  been  made  for  a  Maxwellian 
EEDF  at  an  ionization  fraction  above  10  .  This 

holds  for  the  body  of  the  distribution  but  not  for 
the  high  energy  tail  of  the  EEDF.  While  the 
deviation  shown  in  Fig  9  is  quite  small,  it  is 
nevertheless  very  important  because  the  ionization 
rate  is  controlled  by  those  electrons  with 
energies  above  the  ionization  energy. 


Summarizing  the  points  made  above:  to  describe 
the  kinetics  of  the  low  pressure  discharge,  it  is 
necessary  to  consider  the  spatial  decay  of  the 
cathodic  beam  at  low  currents  and  the  departure 
from  a  Maxwellian  even  at  high  electron  densities. 

When  the  discharge  is  switched  off,  the  ionization 
decays  by  diffusion  and  recombination.  The 
singlet  excited  states  in  H2  decay  down  to,  + 
vibrationally  excited  levels  of  the  H~  (  X  Z.q) 
ground  state.  The  vibrational  level s^appear  to 
have  a  de-excitation  coefficient  of  about  0.1  for 
wall  collisions.  The  electron  temperature  decay 
is  therefore  much  faster  than  the  vibrational 
temperature  decay.  The  electron  energy  equation 
for  a  non-flowing  plasma  after  the  electric  field 
is  removed,  reduces  to 

1L  . tS-Z)  m 


T  is  larger  than  T  however  the  energy 
transfer  coefficient  K  >>  K  ,  therefore  the 
initial  electron  temperatureeaecay  is 

Z-  Tv  +  (ZpT)exp(-^%J)  m 

Therefore  in  some  discharges  after  switch-off, 
the  electron  temperature  will  decay  to  the 
vibrational  temperature  and  the  subsequent  decay 
will  then  depend  on  the  relaxation  time  for 
vibrational  energy.  This  sequence  would  appear  to 
increase  the  diffusion  losses  provided  that 
vibrational  excitation  does  not  cause  enhanced 
dissociative  attachment  in  the  gas  studied. 

The  above  notes  have  treated  in  a  general  manner 
some  non-equilibrium  effects-that  are  important 
when  considering  pulsed  power  switching.  Wave  - 
particle  interaction  and  scattering  will  be 
important  in  some  low  pressure  switches.  The 
electron  kinetics  modeling  assists  in  establishing 
trends  and  limits  of  different  approaches. 

However,  the  most  comprehensive  model  that  does 
not  have  any  assumptions  and  contains  all  of  the 
physical  processes  and  boundary  conditions,  is 
still  a  well-defined  experiment. 
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FIG  1.  Electron  Energy  Distribution  Function 
in  Molecular  Oxygen.  Parameter  is  E/N 
(Townsends) 
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FIG  2.  Electron  Energy  Distribution  Function  in 
Oxygen  50%  Dissociation  (E/N  =  10  Td  and 
8  Td),  0%  Dissociation  (E/N  -  10  Td) 
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Electron  Energy  Distribution  Functions  in 
Hydrogen.  ( E/N  =  40  Td)  Dissociation  as 
Noted. 
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FIG  4.  Electron  Excitation  Rates  in  Molecular 
kV  (01)  First  Vibrational  Level. 
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Influence  of  Vibrational  Temperature  on 
Dissociative  Attachment  Rate  in  Hydrogen. 
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FIG  6.  Calculated  Electron  Drift  Velocity  in 
CO  -  Krypton  Mixtures. 


FIG  7.  Calculated  Electron  Energy  Distribution 
Functions  in  CO  -  Krypton  (2:98) 


FIG  8.  Calculated  Attenuation  and  Scattering  of  Curve)  Dashed  Curve  is  5eV  Maxwellian. 

Fast  Electrons  in  Hydrogen  (0.3  torr)  Modified  Dash:  EQ  =  49.5eV 

Injection  Energy  20eV,  Cathode  Radius  1.5  cm 


